No recombination suppression in asexually produced males of Daphnia pulex.

Obligate parthenogenesis (OP) is often thought to evolve by disruption of reductional meiosis and suppression of crossover recombination. In the crustacean Daphnia pulex, OP lineages, which have evolved from cyclical parthenogenetic (CP) ancestors, occasionally produce males that are capable of reductional meiosis. Here, by constructing high-density linkage maps, we find that these males show only slightly and nonsignificantly reduced recombination rates compared to CP males and females. Both meiosis disruption and recombination suppression are therefore sex-limited (or partly so), which speaks against the evolution of OP by disruption of a gene that is essential for meiosis or recombination in both sexes. The findings may be explained by female-limited action of genes that suppress recombination, but previously identified candidate genes are known to be expressed in both sexes. Alternatively, and equally consistent with the data, OP might have evolved through a reuse of the parthenogenesis pathways already present in CP and through their extension to all events of oogenesis. The causal mutations for the CP to OP transition may therefore include mutations in genes involved in oogenesis regulation and may not necessarily be restricted to genes of the "meiosis toolkit." More generally, our study emphasizes that there are many ways to achieve asexuality, and elucidating the possible mechanisms is key to ultimately identify the genes and traits involved.

Asexual male production by ZW recombination in Artemia parthenogenetica.

In some asexual species, parthenogenetic females occasionally produce males, which may strongly affect the evolution and maintenance of asexuality if they cross with related sexuals and transmit genes causing asexuality to their offspring ("contagious parthenogenesis"). How these males arise in the first place has remained enigmatic, especially in species with sex chromosomes. Here, we test the hypothesis that rare, asexually produced males of the crustacean Artemia parthenogenetica are produced by recombination between the Z and W sex chromosomes during non-clonal parthenogenesis, resulting in ZZ males through loss of heterozygosity at the sex determination locus. We used RAD-sequencing to compare asexual mothers with their male and female offspring. Markers on several sex-chromosome scaffolds indeed lost heterozygosity in all male but no female offspring, suggesting that they correspond to the sex-determining region. Other sex-chromosome scaffolds lost heterozygosity in only a part of the male offspring, consistent with recombination occurring at a variable location. Alternative hypotheses for the production of these males (such as partial or total hemizygosity of the Z) could be excluded. Rare males are thus produced because recombination is not entirely suppressed during parthenogenesis in A. parthenogenetica. This finding may contribute to explaining the maintenance of recombination in these asexuals.

Genetic Drift Shapes the Evolution of a Highly Dynamic Metapopulation.

The dynamics of extinction and (re)colonization in habitat patches are characterizing features of dynamic metapopulations, causing them to evolve differently than large, stable populations. The propagule model, which assumes genetic bottlenecks during colonization, posits that newly founded subpopulations have low genetic diversity and are genetically highly differentiated from each other. Immigration may then increase diversity and decrease differentiation between subpopulations. Thus, older and/or less isolated subpopulations are expected to have higher genetic diversity and less genetic differentiation. We tested this theory using whole-genome pool-sequencing to characterize nucleotide diversity and differentiation in 60 subpopulations of a natural metapopulation of the cyclical parthenogen Daphnia magna. For comparison, we characterized diversity in a single, large, and stable D. magna population. We found reduced (synonymous) genomic diversity, a proxy for effective population size, weak purifying selection, and low rates of adaptive evolution in the metapopulation compared with the large, stable population. These differences suggest that genetic bottlenecks during colonization reduce effective population sizes, which leads to strong genetic drift and reduced selection efficacy in the metapopulation. Consistent with the propagule model, we found lower diversity and increased differentiation in younger and also in more isolated subpopulations. Our study sheds light on the genomic consequences of extinction-(re)colonization dynamics to an unprecedented degree, giving strong support for the propagule model. We demonstrate that the metapopulation evolves differently from a large, stable population and that evolution is largely driven by genetic drift.

Sex ratio effects on reproductive success of male and female Daphnia.

The intensity of mating competition and the opportunity for sexual selection are thought to depend on the operational sex ratio, the ratio of sexually active males to fertilizable females. Cyclic parthenogens, organisms that alternate between sexual reproduction and female-only parthenogenesis, show particularly high variation in sex ratios in natural populations but the effects of this variation on mating competition and reproductive success of each sex are poorly understood. In a series of experiments with Daphnia magna, we experimentally imposed five sex ratio categories, varying from one male per 81 females to an even sex ratio. We found that, in males, reproductive success strongly and monotonically decreased with decreasing number of females per male. In females, in contrast, mating success and reproductive success were reduced only at the most female-biased sex ratio (1:81), when many females remained unmated and unfertilized, and then again at equal sex ratios, probably due to negative effects of high density or stress induced by numerous males. Our results suggest that females experienced male limitation at heavily female-biased sex ratios below one male to about 50 females. As this is well within the sex ratio variation observed in natural Daphnia populations, we conclude that mating competition and the opportunity for sexual selection may exist not only in males but, at least periodically, also in females.

Not so clonal asexuals: Unraveling the secret sex life of Artemia parthenogenetica.

The maintenance of sex is paradoxical as sexual species pay the "twofold cost of males" and should thus quickly be replaced by asexual mutants reproducing clonally. However, asexuals may not be strictly clonal and engage in "cryptic sex," challenging this simple scenario. We study the cryptic sex life of the brine shrimp Artemia parthenogenetica, which has once been termed an "ancient asexual" and where no genetic differences have ever been observed between parents and offspring. This asexual species rarely produces males, which can hybridize with sexual females of closely related species and transmit asexuality to their offspring. Using such hybrids, we show that recombination occurs in asexual lineages, causing loss-of-heterozygosity and parent-offspring differences. These differences cannot generally be observed in field-sampled asexuals because once heterozygosity is lost, subsequent recombination leaves no footprint. Furthermore, using extensive paternity tests, we show that hybrid females can reproduce both sexually and asexually, and transmit asexuality to both sexually and asexually produced offspring in a dominant fashion. Finally, we show that, contrary to previous reports, field-sampled asexual females also rarely reproduce sexually (rate ∼2‰). Overall, most previously known facts about Artemia asexuality turned out to be erroneous. More generally, our findings suggest that the evidence for strictly clonal reproduction of asexual species needs to be reconsidered, and that rare sex and consequences of nonclonal asexuality, such as gene flow within asexuals, need to be more widely taken into account in more realistic models for the maintenance of sex and the persistence of asexual lineages.

Parasite-driven replacement of a sexual by a closely related asexual taxon in nature.

Asexual species are thought to suffer more from coevolving parasites than related sexuals. Yet a variety of studies do not find the patterns predicted by theory. Here, to shine light on this conundrum, we investigate one such case of an asexual advantage in the presence of parasites. We follow the frequency dynamics of sexual and asexual Daphnia pulex in a natural pond that was initially dominated by sexuals. Coinciding with an epidemic of a microsporidian parasite infecting both sexuals and asexuals, the pond was rapidly taken over by the initially rare asexuals. With experiments comparing multiple sexual and asexual clones from across the local metapopulation, we confirm that asexuals are less susceptible and also suffer less from the parasite once infected. These results are consistent with the parasite-driven, ecological replacement of dominant sexuals by closely related, but more resistant asexuals, ultimately leading to the extinction of the formerly superior sexual competitor. Our study is one of the clearest examples from nature, backed up by experimental verification, showing a parasite-mediated reversal of competition dynamics. The experiments show that, across the metapopulation, asexuals have an advantage in the presence of parasites. In this metapopulation, asexuals are relatively rare, likely due to their recent invasion. While we cannot rule out other reasons for the observed patterns, the results are consistent with a temporary parasite-mediated advantage of asexuals due to the fact that they are rare, which is an underappreciated aspect of the Red Queen Hypothesis.

Genetic control of male production in Daphnia pulex.

Daphnia normally reproduce by cyclical parthenogenesis, with offspring sex being determined by environmental cues. However, some females have lost the ability to produce males. Our results demonstrate that this loss of male-producing ability is controlled by a dominant allele at a single locus. We identified the locus by comparing whole-genome sequences of 67 nonmale-producing (NMP) and 100 male-producing (MP) clones from 5 Daphnia pulex populations, revealing 132 NMP-linked SNPs and 59 NMP-linked indels within a single 1.1-Mb nonrecombining region on chromosome I. These markers include 7 nonsynonymous mutations, all of which are located within one unannotated protein-coding gene (gene 8960). Within this single gene, all of the marker-linked NMP haplotypes from different populations form a monophyletic clade, suggesting a single origin of the NMP phenotype, with the NMP haplotype originating by introgression from a sister species, Daphnia pulicaria Methyl farnesoate (MF) is the innate juvenile hormone in daphnids, which induces the production of males and whose inhibition results in female-only production. Gene 8960 is sensitive to treatment by MF in MP clones, but such responsiveness is greatly reduced in NMP clones. Thus, we hypothesize that gene 8960 is located downstream of the MF-signaling pathway in D. pulex, with the NMP phenotype being caused by expression change of gene 8960.

Evolution of Gene Expression during a Transition from Environmental to Genetic Sex Determination.

Genetic sex determination (GSD) can evolve from environmental sex determination (ESD) via an intermediate state in which both coexist in the same population. Such mixed populations are found in the crustacean Daphnia magna, where non-male producers (NMP, genetically determined females) coexist with male producers (MP), in which male production is environmentally inducible and can also artificially be triggered by exposure to juvenile hormone. This makes Daphnia magna a rare model species for the study of evolutionary transitions from ESD to GSD. Although the chromosomal location of the NMP-determining mutation has been mapped, the actual genes and pathways involved in the evolution of GSD from ESD remain unknown. Here, we present a transcriptomic analysis of MP and NMP females under control (female producing) and under hormone exposure conditions. We found ∼100 differentially expressed genes between MP and NMP under control conditions. Genes in the NMP-determining chromosome region were especially likely to show such constitutive expression differences. Hormone exposure led to expression changes of an additional ∼100 (MP) to ∼600 (NMP) genes, with an almost systematic upregulation of those genes in NMP. These observations suggest that the NMP phenotype is not determined by a simple "loss-of-function" mutation. Rather, homeostasis of female offspring production under hormone exposure appears to be an active state, tightly regulated by complex mechanisms involving many genes. In a broader view, this illustrates that the evolution of GSD, while potentially initiated by a single mutation, likely leads to secondary integration involving many genes and pathways.

How clonal are clones? A quest for loss of heterozygosity during asexual reproduction in Daphnia magna.

Due to the lack of recombination, asexual organisms are predicted to accumulate mutations and show high levels of within-individual allelic divergence (heterozygosity); however, empirical evidence for this prediction is largely missing. Instead, evidence of genome homogenization during asexual reproduction is accumulating. Ameiotic crossover recombination is a mechanism that could lead to long genomic stretches of loss of heterozygosity (LOH) and unmasking of mutations that have little or no effect in heterozygous state. Therefore, LOH might be an important force for inducing variation among asexual offspring and may contribute to the limited longevity of asexual lineages. To investigate the genetic consequences of asexuality, here we used high-throughput sequencing of Daphnia magna for assessing the rate of LOH over a single generation of asexual reproduction. Comparing parthenogenetic daughters with their mothers at several thousand genetic markers generated by restriction site-associated DNA (RAD) sequencing resulted in high LOH rate estimation that largely overlapped with our estimates for the error rate. To distinguish these two, we Sanger re-sequenced the top 17 candidate RAD-loci for LOH, and all of them proved to be false positives. Hence, even though we cannot exclude the possibility that short stretches of LOH occur in genomic regions not covered by our markers, we conclude that LOH does not occur frequently during asexual reproduction in D. magna and ameiotic crossovers are very rare or absent. This finding suggests that clonal lineages of D. magna will remain genetically homogeneous at least over time periods typically relevant for experimental work.

Phylogeography of Daphnia magna Straus (Crustacea: Cladocera) in Northern Eurasia: Evidence for a deep longitudinal split between mitochondrial lineages.

Species with a large geographic distributions present a challenge for phylogeographic studies due to logistic difficulties of obtaining adequate sampling. For instance, in most species with a Holarctic distribution, the majority of studies has concentrated on the European or North American part of the distribution, with the Eastern Palearctic region being notably understudied. Here, we study the phylogeography of the freshwater cladoceran Daphnia magna Straus, 1820 (Crustacea: Cladocera), based on partial mitochondrial COI sequences and using specimens from populations spread longitudinally from westernmost Europe to easternmost Asia, with many samples from previously strongly understudied regions in Siberia and Eastern Asia. The results confirm the previously suspected deep split between Eastern and Western mitochondrial haplotype super-clades. We find a narrow contact zone between these two super-clades in the eastern part of Western Siberia, with proven co-occurrence in a single lake in the Novosibirsk region. However, at present there is no evidence suggesting that the two mitochondrial super-clades represent cryptic species. Rather, they may be explained by secondary contact after expansion from different refugia. Interestingly, Central Siberia has previously been found to be an important contact zone also in other cladoceran species, and may thus be a crucial area for understanding the Eurasian phylogeography of freshwater invertebrates. Together, our study provides an unprecedented complete, while still not global, picture of the phylogeography of this important model species.

Identification of General Patterns of Sex-Biased Expression in Daphnia, a Genus with Environmental Sex Determination.

Daphnia reproduce by cyclic-parthenogenesis, where phases of asexual reproduction are intermitted by sexual production of diapause stages. This life cycle, together with environmental sex determination, allow the comparison of gene expression between genetically identical males and females. We investigated gene expression differences between males and females in four genotypes of Daphnia magna and compared the results with published data on sex-biased gene expression in two other Daphnia species, each representing one of the major phylogenetic clades within the genus. We found that 42% of all annotated genes showed sex-biased expression in D. magna This proportion is similar both to estimates from other Daphnia species as well as from species with genetic sex determination, suggesting that sex-biased expression is not reduced under environmental sex determination. Among 7453 single copy, one-to-one orthologs in the three Daphnia species, 707 consistently showed sex-biased expression and 675 were biased in the same direction in all three species. Hence these genes represent a core-set of genes with consistent sex-differential expression in the genus. A functional analysis identified that several of them are involved in known sex determination pathways. Moreover, 75% were overexpressed in females rather than males, a pattern that appears to be a general feature of sex-biased gene expression in Daphnia.

Low recombination rates in sexual species and sex-asex transitions.

In most sexual, diploid eukaryotes, at least one crossover occurs between each pair of homologous chromosomes during meiosis, presumably in order to ensure proper segregation. Well-known exceptions to this rule are species in which one sex does not recombine and specific chromosomes lacking crossover. We review other possible exceptions, including species with chromosome maps of less than 50 cM in one or both sexes. We discuss the idea that low recombination rates may favour sex-asex transitions, or, alternatively may be a consequence of it. We then show that a yet undescribed species of brine shrimp Artemia from Kazakhstan (A sp. Kazakhstan), the closest known relative of the asexual Artemia parthenogenetica, has one of the shortest genetic linkage maps known. Based on a family of 42 individuals and 589 RAD markers, we find that many linkage groups are considerably shorter than 50 cM, suggesting either no obligate crossover or crossovers concentrated at terminal positions with little effect on recombination. We contrast these findings with the published map of the more distantly related sexual congener, A. franciscana, and conclude that the study of recombination in non-model systems is important to understand the evolutionary causes and consequences of recombination.This article is part of the themed issue 'Evolutionary causes and consequences of recombination rate variation in sexual organisms'.

Transition from Environmental to Partial Genetic Sex Determination in Daphnia through the Evolution of a Female-Determining Incipient W Chromosome.

Sex chromosomes can evolve during the evolution of genetic sex determination (GSD) from environmental sex determination (ESD). Despite theoretical attention, early mechanisms involved in the transition from ESD to GSD have yet to be studied in nature. No mixed ESD-GSD animal species have been reported, except for some species of Daphnia, small freshwater crustaceans in which sex is usually determined solely by the environment, but in which a dominant female sex-determining locus is present in some populations. This locus follows Mendelian single-locus inheritance, but has otherwise not been characterized genetically. We now show that the sex-determining genomic region maps to the same low-recombining peri-centromeric region of linkage group 3 (LG3) in three highly divergent populations of D. magna, and spans 3.6 Mb. Despite low levels of recombination, the associated region contains signs of historical recombination, suggesting a role for selection acting on several genes thereby maintaining linkage disequilibrium among the 36 associated SNPs. The region carries numerous genes involved in sex differentiation in other taxa, including transformer2 and sox9. Taken together, the region determining the genetic females shows characteristics of a sex-related supergene, suggesting that LG3 is potentially an incipient W chromosome despite the lack of significant additional restriction of recombination between Z and W. The occurrence of the female-determining locus in a pre-existing low recombining region illustrates one possible form of recombination suppression in sex chromosomes. D. magna is a promising model for studying the evolutionary transitions from ESD to GSD and early sex chromosome evolution.

A high-density genetic map reveals variation in recombination rate across the genome of Daphnia magna.

BACKGROUND: Recombination rate is an essential parameter for many genetic analyses. Recombination rates are highly variable across species, populations, individuals and different genomic regions. Due to the profound influence that recombination can have on intraspecific diversity and interspecific divergence, characterization of recombination rate variation emerges as a key resource for population genomic studies and emphasises the importance of high-density genetic maps as tools for studying genome biology. Here we present such a high-density genetic map for Daphnia magna, and analyse patterns of recombination rate across the genome. RESULTS: A F2 intercross panel was genotyped by Restriction-site Associated DNA sequencing to construct the third-generation linkage map of D. magna. The resulting high-density map included 4037 markers covering 813 scaffolds and contigs that sum up to 77 % of the currently available genome draft sequence (v2.4) and 55 % of the estimated genome size (238 Mb). Total genetic length of the map presented here is 1614.5 cM and the genome-wide recombination rate is estimated to 6.78 cM/Mb. Merging genetic and physical information we consistently found that recombination rate estimates are high towards the peripheral parts of the chromosomes, while chromosome centres, harbouring centromeres in D. magna, show very low recombination rate estimates. CONCLUSIONS: Due to its high-density, the third-generation linkage map for D. magna can be coupled with the draft genome assembly, providing an essential tool for genome investigation in this model organism. Thus, our linkage map can be used for the on-going improvements of the genome assembly, but more importantly, it has enabled us to characterize variation in recombination rate across the genome of D. magna for the first time. These new insights can provide a valuable assistance in future studies of the genome evolution, mapping of quantitative traits and population genetic studies.

Evolutionary mysteries in meiosis.

Meiosis is a key event of sexual life cycles in eukaryotes. Its mechanistic details have been uncovered in several model organisms, and most of its essential features have received various and often contradictory evolutionary interpretations. In this perspective, we present an overview of these often 'weird' features. We discuss the origin of meiosis (origin of ploidy reduction and recombination, two-step meiosis), its secondary modifications (in polyploids or asexuals, inverted meiosis), its importance in punctuating life cycles (meiotic arrests, epigenetic resetting, meiotic asymmetry, meiotic fairness) and features associated with recombination (disjunction constraints, heterochiasmy, crossover interference and hotspots). We present the various evolutionary scenarios and selective pressures that have been proposed to account for these features, and we highlight that their evolutionary significance often remains largely mysterious. Resolving these mysteries will likely provide decisive steps towards understanding why sex and recombination are found in the majority of eukaryotes.This article is part of the themed issue 'Weird sex: the underappreciated diversity of sexual reproduction'.

Daphnia magna transcriptome by RNA-Seq across 12 environmental stressors.

The full exploration of gene-environment interactions requires model organisms with well-characterized ecological interactions in their natural environment, manipulability in the laboratory and genomic tools. The waterflea Daphnia magna is an established ecological and toxicological model species, central to the food webs of freshwater lentic habitats and sentinel for water quality. Its tractability and cyclic parthenogenetic life-cycle are ideal to investigate links between genes and the environment. Capitalizing on this unique model system, the STRESSFLEA consortium generated a comprehensive RNA-Seq data set by exposing two inbred genotypes of D. magna and a recombinant cross of these genotypes to a range of environmental perturbations. Gene models were constructed from the transcriptome data and mapped onto the draft genome of D. magna using EvidentialGene. The transcriptome data generated here, together with the available draft genome sequence of D. magna and a high-density genetic map will be a key asset for future investigations in environmental genomics.

Genetic load, inbreeding depression, and hybrid vigor covary with population size: An empirical evaluation of theoretical predictions.

Reduced population size is thought to have strong consequences for evolutionary processes as it enhances the strength of genetic drift. In its interaction with selection, this is predicted to increase the genetic load, reduce inbreeding depression, and increase hybrid vigor, and in turn affect phenotypic evolution. Several of these predictions have been tested, but comprehensive studies controlling for confounding factors are scarce. Here, we show that populations of Daphnia magna, which vary strongly in genetic diversity, also differ in genetic load, inbreeding depression, and hybrid vigor in a way that strongly supports theoretical predictions. Inbreeding depression is positively correlated with genetic diversity (a proxy for Ne ), and genetic load and hybrid vigor are negatively correlated with genetic diversity. These patterns remain significant after accounting for potential confounding factors and indicate that, in small populations, a large proportion of the segregation load is converted into fixed load. Overall, the results suggest that the nature of genetic variation for fitness-related traits differs strongly between large and small populations. This has large consequences for evolutionary processes in natural populations, such as selection on dispersal, breeding systems, ageing, and local adaptation.

Uncovering Cryptic Asexuality in Daphnia magna by RAD Sequencing.

The breeding systems of many organisms are cryptic and difficult to investigate with observational data, yet they have profound effects on a species' ecology, evolution, and genome organization. Genomic approaches offer a novel, indirect way to investigate breeding systems, specifically by studying the transmission of genetic information from parents to offspring. Here we exemplify this method through an assessment of self-fertilization vs. automictic parthenogenesis in Daphnia magna. Self-fertilization reduces heterozygosity by 50% compared to the parents, but under automixis, whereby two haploid products from a single meiosis fuse, the expected heterozygosity reduction depends on whether the two meiotic products are separated during meiosis I or II (i.e., central vs. terminal fusion). Reviewing the existing literature and incorporating recombination interference, we derive an interchromosomal and an intrachromosomal prediction of how to distinguish various forms of automixis from self-fertilization using offspring heterozygosity data. We then test these predictions using RAD-sequencing data on presumed automictic diapause offspring of so-called nonmale producing strains and compare them with "self-fertilized" offspring produced by within-clone mating. The results unequivocally show that these offspring were produced by automixis, mostly, but not exclusively, through terminal fusion. However, the results also show that this conclusion was only possible owing to genome-wide heterozygosity data, with phenotypic data as well as data from microsatellite markers yielding inconclusive or even misleading results. Our study thus demonstrates how to use the power of genomic approaches for elucidating breeding systems, and it provides the first demonstration of automictic parthenogenesis in Daphnia.

Genes mirror geography in Daphnia magna.

Identifying the presence and magnitude of population genetic structure remains a major consideration in evolutionary biology as doing so allows one to understand the demographic history of a species as well as make predictions of how the evolutionary process will proceed. Next-generation sequencing methods allow us to reconsider previous ideas and conclusions concerning the distribution of genetic variation, and what this distribution implies about a given species evolutionary history. A previous phylogeographic study of the crustacean Daphnia magna suggested that, despite strong genetic differentiation among populations at a local scale, the species shows only moderate genetic structure across its European range, with a spatially patchy occurrence of individual lineages. We apply RAD sequencing to a sample of D. magna collected across a wide swath of the species' Eurasian range and analyse the data using principle component analysis (PCA) of genetic variation and Procrustes analytical approaches, to quantify spatial genetic structure. We find remarkable consistency between the first two PCA axes and the geographic coordinates of individual sampling points, suggesting that, on a continent-wide scale, genetic differentiation is driven to a large extent by geographic distance. The observed pattern is consistent with unimpeded (i.e. no barriers, landscape or otherwise) migration at large spatial scales, despite the fragmented and patchy nature of favourable habitats at local scales. With high-resolution genetic data similar patterns may be uncovered for other species with wide geographic distributions, allowing an increased understanding of how genetic drift and selection have shaped their evolutionary history.